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Abstract 

We calculate the magnetization dynamics induced by the inverse Faraday effect in disordered 
metals in THz regime by using the diagrammatic method. We find that the induced magnetization 
is proportional to the frequency of circularly polarized light. 



I. INTRODUCTION 



Photoinduced magnetization by circularly polarized light (the inverse Faraday effect) has 
been studied since its theoretical prediction in the 1960s- - -. Recently, experiments using 
intense laser have been conducted, and the induced effective magnetic field was shown to be 
as strong as a few TeslaA~-. The fast magnetization reversal realized is expected to be useful 
for application to magnetic devices^. 

The mechanism of the inverse Faraday effect was discussed by Pitaevskii based on a 
symmetry argument^. He noted that the free energy of the electromagnetic field in solids is 
generally given by 



where is a 3 x 3 tensor of the dielectric functions dependent on mediums, and Ei is the 
complex representation of the electric field. The magnetization induced by light is given by 
the functional derivative of F with respect to the external magnetic field, H ex , as 



Because of the Onsager relation, the asymmetric components of are linear or higher odd 
power in the external magnetic field, while the symmetric components are an even power 
of the magnetic fieldi^. Therefore, expanding up to the linear order in the magnetic field, 
the dielectric tensor is = e^Sij + ixi^ijkH^ x , where Eq and xi are constants independent 
of i? ex . The induced magnetization is therefore written generally as 



It is thus proportional to the helicity vector of the circularly polarized light, i(E x E*). 
(E x E* is parallel to light propagation and is pure imaginary) 1-5 . Because the helicity 
vector, i{E x E*), of the oscillate phase factor is canceled out for Ee iQt and E*e~ tQt , M 
is a static and rectified magnetization. This magnetization is converted from the optical 
angular momentum of circularly polarized light, E x E*, and this conversion efficiency 
is determined by xi- The coefficient xi was explicitly calculated in the visible light case 
considering plasma oscillation^. It was found to strongly depend on the angular frequency Q 
as Xi °c f2~ 3 . Another mechanism of the inverse Faraday effect is spin-orbit interaction, as 




(1) 





M = %X\E x E*. 



(3) 
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discussed by Pershan and coworkers^-. They also showed that the spin-orbit contribution 
also results in xi ^ ^ 3 - 

Besides the effect of the visible light lasers, the effect of terahertz lasers on magnetic 
systems has recently been attracting interest. For example, the spin structure and excitation 
in ferroelectric magnets were explored by utilizing use of the magnetoelectric effect due to 
a terahertz laser—. 

In this paper, we will theoretically investigate the magnetization induced by the inverse 
Faraday effect in a THz regime. The magnetization is calculated using the diagrammatic 
method in disordered metals considering the spin-orbit interaction. The induced magneti- 
zation is found to linearly depend on the frequency in the THz regime. 

II. HAMILTON! AN 

We consider the conduction electrons in the presence of spin-orbit interaction and im- 
purity scattering. The applied THz light is represented by the vector potential A em . The 
Hamiltonian of the system is % = Hq + T-L em + ?4o + Hi, where 



and c, a are the conduction electron's annihilation and creation operators, respectively. We 
consider two origins for spin-orbit interaction: one from the applied electromagnetic field 
and the other induced by random impurities. The coefficient A = 7-^-3 is the strength of 
spin-orbit coupling, and U{ — u J2e$( x ~ ^) ^ s ^ ne impurity potential, where Ri represents 
the randomly distribution impurity positions. (The average of Am represents the strength of 
the spin-orbit interaction in the material.) The scattering of the electron by the impurities 
is represented by We consider the elastic lifetime of the electron, r, arising from this 
scattering. 

We consider a monochromatic light with a frequency Q and wave vector q(The unit vector 
of electromagnetic field, q = q/\q\, is parallel to E x E*). The vector potential is related 





(7) 



(4) 
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to the electric field by £ cm = —A cm . We express the vector potential and electric field by 
utilizing the complex amplitudes E and A as 

£ cm = l(Ee i{q - x - nt) + E*e- i{q - x ' nt) ), (8) 
A cm = ]-(Ae i( - q - x - nt) + A*e- i ^ x - nt) ). (9) 

The amplitude thus satisfies A = — and A* = i^-. 



III. CALCULATION OF THE MAGNETIZATION DUE TO THE INVERSE 
FARADAY EFFECT 

The spin density induced by the applied circularly polarized light (inverse Faraday ef- 
fect) is calculated by estimating the expectation value s a = (cV a c) (a = x,y,z represents 
the spin direction, and ( ) denotes expectation values). Using the lesser Green's function, 
G < (x,x;t,t), the spin density is expressed 

s a (x,t) = -ihTr[a a G < (x,x;t,t)}. (10) 

The calculation of the spin density is performed by treating "H cm and "H so perturbatively. 
From the spin density, the magnetization of inverse Faraday effect is given by M a = —^s a , 
where /zb is the Bohr magneton, and g is the g-factor. In the following section, we only con- 
sider the contribution proportional to Ex E*, because we are interested in the magnetization 
induced by the inverse Faraday effect. 



A. Spin-orbit interaction from the electric field 

We first estimate the spin-orbit interaction arising from the applied electric field, whose 
Hamiltonian we represent as "Hso\ i.e., 

H£ } =~\j dx3c ^ ■ Mem X (p - eA cm )]c. (11) 

The contribution to the spin density, s^ a , due to l-i^o are shown in Fig. 1, and are calculated 

as 

s (i)« = - t 2e 2 \e aye {E 1 A\ + E*A t ) K(Sl), (12) 
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FIG. 1: Diagrams contributing to the spin density sW induced by circularly polarized light. The 
electron's Green's functions including the impurity average here are denoted by thick lines. The 
wavy line and the dashed arrow represent the gauge field of light (A em ) and electric field (A em ), 
respectively. A is spin-orbit coupling constant accompanied with A cm . The vertex shows the Pauli 
matrix (<t). 



where 



K{0) = 

k,ui L 

m 



H 



in 



(13) 



= — S^k t k c [g 

nm ' ' 



fc,U) 9k,U}+fl 9k,ui\ 



k,ui 



and gk tLJ is the free Green's function on a Keldysh contour. The lesser component of the free 
Green's function is g£ u = f{oj){g% w — 9kJ)i where f(cu) is the Fermi distribution function, 
g a = (hw — €k — i^-)" 1 and g r = (g a )* are the advanced and retarded Green's functions, 
respectively^, (e^ = ^p- — and eF is the Fermi energy). Here K(£l) is the even function 
with respect to Q. 

In the terahertz regime, the frequency of light, Q, is much smaller than that of the 
electron, u>. We therefore expand the Fermi distribution function as f(u + Q) = f(u) + 
Qf(u) + -¥f"(u)) + o{Q 3 ). We thus obtain 

k,u) 

where we have used the partial integration with respect to k (Note that i^(0) = 0). Con- 
sidering zero temperature, i.e., f'ioj) = —5(u), we obtain K(Q) as 



K(tt) 



244 



hn 2 5 ec + o(n 4 ). 



From Eqs. ( fl~2l) and ( JT4l) . the result of the induced spin density is given by 

7r^e 2 A 



6e F 



-hQeap^EpE*, 



(14) 



(15) 
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where v is the electron density of state at the Fermi energy per volume. 



B. Spin-orbit interaction due to impurities 

Next, we consider the inverse Faraday effect from spin-orbit interaction caused by the 
random impurity potential: 

Hg = ^J dx 3 Ja ■ [V«i x (p - eA cm )]c. (16) 

The random potential is treated by averaging, {ui(q)ui(q')) i = n{uf5(q + q') and 
(u\(q)u\(q')ui(q"))i = riiuf5(q + q' + q"), where q, q', q" are wave numbers, ri\ is the concen- 
tration of the impurity, and ( )j is the random impurity average. 

Contributions to the spin density proportional to E x E* are diagrammatically shown 
in Figs. |2] and [3j Processes shown in Fig. [2] and Fig. [3] differ in impurity averaging; the 
averaging in Fig. [2] is third order in the impurities as in the case of the anomalous Hall 
effect^. We first estimate the contribution (a) of Fig. El defined as sf^. It is 

s 3(a) = ^—e^EeE^Ta - r_ n ). (17) 

In the above equation, we have used the relation, A cm = | (^ e ' tQt — ^ e%m )- is written 
as 



■p 2 \ -* 2 , . / 9k,cu9k,u! (.9k,u)+n9k',u)+n 9k",ui+n C - C J \ 

k,k',k",cu \~9k,uj+n 9k,ui+n\9k,ui 9w,w 9h",w ~ c - c 



n 



2 

k,k'k",w 



2 r / \ 



(18) 



and £r is given by 



2 ^ 

^ n = ?)2 ^ 2 ^ aJ (fl'fc'.w 9k" ) u{9k,u + fl'fe'.w + ^fc",w) — 9k,ui9k',uj+n 9k",w+n)) ~ c - c ] 



n 2 

k,k'k",UI 



(19) 

The first term of Tq is even function respect with Q, and therefore Tq — T_q = 5T^ — 5T_q. 
The Green's function can expanded as g a k u)+n = gl Ld -Ml(glJ 2 + (m) 2 (glJ 2 -(hn) 3 (glJ 3 + 
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FIG. 2: Diagrams contributing to the spin density in the spin-orbit interaction due to random 
impurities. Dashed lines denote the impurity averaging between the spin-orbit interaction(the 
combination between the wavy line and the open circle) and the impurity scattering (open circle). 
The contribution of diagrams (a) and (b) indicates sj^? and si^f, respectively. 

o(f2 4 ), and thus 5Tn — ST^q is reduced to 



This contribution at low frequency contains only the terms of considering of only the retarded 
or advanced Green's functions, and it turns out to be negligibly small compared with Jq of 



5T n -5T_ n = 4h 3 n ^ 



k,k',k",u 



( +(^,J 3 (^,J 4 <fc + (gV 3 9 a k>,u (g a k »J 4 \ 

+K J 4 te,J 3 <fc + (g a k> JYk>,„ J 3 
+(gtJ 3 (g a k >,J 2 J 3 + (#J 8 te2>) 8 (g a k »,J 2 
V +(gt,J\g a k >J 2 (g a k »,J 2 J 



— c.c 



+ o(n 3 ). 



(20) 



Eq. ([22]). 

The contribution of Fig. 12(b), s. 



,(2)a 
3(b) ' 



IS 
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where 



Jc 



h 2 



3mn 2 



FIG. 3: Diagrammatic representation of spin density. 



^ (k') 2 (k") 2 [gk tU > gk",uj <7fc",u>+n gk',ui+n gk',u 5'fc,w] < ^ 
+k 2 {k") 2 [gk,ui gk,ui+n gk»,u)+n gk",w gk',ui <?fc,w] 

+k 2 (k») 2 [g k ,u) gk",u) gk",u)+a gk',Lu+a gk,w+si gk,w\ 

+k 2 {k') 2 [g k 

,u} g k ",Lu g k ' ,ujg k ' ,u>+£i g k ,uj+Q gk,u\ 
\+k 2 (k') 2 [g k g k ,w+n gk",uj+n g k ',uj+n g k ',ui g k ,J\ J 



k,k'k"w 



The dominant contribution of Jn is estimated as 

/ 



Jo 



h 2 



£ k 2 {k"f{Un-f.) 



k,k\k'\ui 



\ 



V 



Wk',J\ 9k,Log k ",u}(9 k ^+n9k",uj+n ~ c - c ) 

+ Iffc,^! 2 (g k ' ,w9k" ,oj + 9l' ,uj9k" ,w)(gk,uj+ngk" — c -°) 
■\ 2 (g k ",Lu + gl" ,Lu)(gk,iu+Qgk' ,u)+ng k " ,uj+n ~ c -°) j 

(22) 



+ 



where 



/ 



SJc, 



3mQ 2 



k,k',k".u) 



\ 



C.C. 



(23) 



(gt' ,lj) 2 9 k ,u}9k" ^gk^+ngt" ,lu+q 

+2 (gk jU ) 2 gki,u9k" ,w9 k ,Lu+v,9k",Lu+n 

^(gk,u) 2 g^ ,^gk,uj+ng k ' ,uj+v,gk'' ,uj+n J 

Because 5Tq and 5Jn have only the contribution only retarded or advanced Green's func- 
tions, Eqs. ( HI?]) and (|2"3"|) are both negligibly small compared with Eq. f[2"2l (by the order of 



<C 1, where r = 



2nvnu 2 



is a elastic electron lifetime). From these estimations, the spin 
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density contribution shown in Fig. El sjj a = sf^ + sf^, is finally obtained as 



S ( 2 ) Q —i3e 2 h\riiU 3 e a e 1 QEgE 1 



r|2 
9k\ 



= il27r 2 eVA (—- J Qre^E^E* 

The contribution of Fig. |3] is calculated similarly as 

(2)q e 2 h 2 \n\u 2 



(24) 



3m 



-Cae-yEpE^ [Iq — J_ 



(25) 



where 



-fo = ^ (^') 2 [S^w 9k', w 9k',io+n 9k,w + 9k,ui 9k',w+n 9k', w 9k,i 



k,k'u) 



+ ^ 2 ^ fc ' w ^ fe/ > w 9k,uj+n 9k,u + S'fc.w 9k,uj+n 9k 1 , u fl , fc,w] < 

^02 ^ 2 [^ fc > w 9k',ui+n 9k,u+n 9k,uj + ^fc.w 9k,w+n 9k>,ui+n 9k,J\ 



k,k'u) 



The odd function of Q, In — I-a, is calculated as 

/ 



2 x - 

In - I-n = 1^ k 



r 9k r 9i 

,o)+f2 \9k,u> 9k',u> 9k,u> 9k',oj) 

l ~(fui+n — fu)9k,w 9k,u (.9k,u,+n 9k',u+ti ~ 9k,u+n 9k',w+n) t 
Ah 3 Q H,J 3 (9t',J 3 GC + 9 a k ',J - c.c] + o(fi 4 ) 



k,k'u) 



(26) 



k,k'u> 



2e* 



+ o(0 J 



(27) 



From Eqs. (125]) and (1271) . reads 



,(2)a 
6 2 



Z7T 2 e 2 \h>ftTL ( h x 



12 4 



7' 



(28) 



Equation (1281) is small compared with Eq. (I24p by the order of ^ (efr) ^ ^ anc ^ thus 

(2) (2) 

neglected below. Therefore, the dominant contribution of w so is given by S3 . 

We summarize the spin density contributed from the spin-orbit interaction and 
'HP'}. Since are negligible small compared with from Eq. ( |T5l) and ( 1281) . the total 
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magnetization of the inverse Faraday effect in spin-orbit interaction, M ~ —^f-s\ is 
obtained as 

M ~ % X E x E*, (29) 

where 

X = -i6gfi B n 2 e 2 v 2 \ (^) fir (30) 

The magnetization is proportional to the strength of spin-orbit coupling, intensity of light, 
frequency, and electron's relaxation time. 

IV. SUMMARY AND DISCUSSION 

We have shown that the spin density induced by the inverse Faraday effect in the THz 
regime is proportional to E x E*. From the spin density induced by circularly polarized 
light, the magnitude of the magnetic field generated in the medium, B eS can be estimated as 
B cS = [iM, where \x are magnetic susceptibility in the medium. Here, we chose the magnetic 
susceptibility in the vacuum. The magnitude of spin-orbit coupling is \k\ = 0.01 — 1 for 
various materials^, and we choose here Xkp ~ 0.5 considering the case of platinum^. We 
choose the amplitude and frequency of applied electromagnetic field as \E\ = 10 8 V/m and 
f2 = ITHz, respectively. In metals with eF ~ 1 eV, ^ = 0.1, and ^ ~ 0.01, the magnetic 
field can be estimated by \B cS \ f» 10~ 2 T. We have theoretically studied the spin density 
induced by circularly polarized light in THz regime in metals spin-orbit interaction. The 
induced spin is proportional to spin-orbit interaction, frequency of applied THz light, and 
also depends on the electrons relaxation time. 
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